Possible fractal structure of exact universal exchange-correlation potential 
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Abstract: Heuristic arguments are presented to show that the basic quantity of density func- 
tional theory - exact universal exchange-correlation potential should have fractal structure. 



The fundamental equation describing electronic structure of matter is N-electron Schrodinger's equation 



T + V + U 
where the kinetic energy operator 
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electron-nuclei interaction potential 
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the sum on k run over all nuclei in the system, each with charge Z k e and position And for a Coulomb system we 
have electron-electron interaction potential 
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The usual quantum-mechanical approach to Schrodinger's equation (SE) can be summarized by the following sequence 

v(r) =^I> ^(ri, r 2 . . . , rjv) observables, (5) 

i.e., one specifies the system by choosing v(r), put it into Schrodinger's equation, solves that equation for the wave 
function ty, and then calculates expectation values of observables with this wave function. Among the observables 
that are calculated in this way is also density 



n(r)=N J d 3 r 2 J d 3 r 3 . . . J d 3 r N ^* (r,r 2 ... ,r N )^(r,r 2 ... ,r N ). 
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Many sophisticated methods for solving Schrodinger's equation have been developed, for example, diagrammatic 
perturbation theory [1] , coupled cluster method [2] , or Brillouin-Wigner coupled cluster method [3] . 

Nevertheless there is an attractive alternative to these methods, the density-functional approach (DFT) [4-7]. 

At the heart of DFT is the Hohenberg-Kohn theorem [9] . This theorem states that for ground states Eq. (6) can 
be inverted: given a ground-state density no(r) it is possible, in principle, to calculate the corresponding ground-state 
wave function v E , o( r i I r 2 ■ ■ ■ , *n)- This means that is a functional of uq. 

Recalling that the kinetic and interaction energies of a nonrelativistic Coulomb system are described by universal 
operators, we can N-electron energy functional E v write as 



E v [n] =T[n] + U[n] + V[n], 
where T[n] and U[n] are universal Junctionals. On the other hand, 



V[n] 
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is obviously nonuniversal (it depends on v(r), i.e., on the system under study), but simple: once the system is specified, 
i.e., v(r) is known, the functional V[n] is known explicitly. All one has to do is to minimize £7„[n] with respect to 
n(r). The minimizing function rio(r) is the system's ground state charge density and the value E Vt o = E v [no] is the 
ground state energy The problem is that the minimization of E v [n] is that the exact functional is not known and one 
needs to used approximations for T[n] and U[n\. 

In the Kohn-Sham approach the kinetic energy T[n], is expressed in terms of the single-particle orbitals <f>i(r) of a 
noninteracting system with density n, as 

h 2 N r 

T *N = "2^E / ^rtf (r)V 2 &(r), (9) 
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because for noninteracting particles the total kinetic energy is just the sum of the individual kinetic energies. Since 
all 0i(r) are functional of n, this expression for T s is an explicit orbital functional but an implicit density functional, 
T s [n}=T s [{<Pi[n}}}. 

We now rewrite the exact energy functional as 

E[n] = T s [n] + U[n] + V[n] = T s [{^[n}}} + U H [n] + E xc [n] + V[n], (10) 

where E xc [n], is called the exchange-correlation (xc) energy functional, which incorporates correlations due to the 
Pauli exclusion principle, Coulomb repulsion, and also correlation contribution to the kinetic energy not included in 
Eq. (9). 

Since T s is now written as an orbital functional one cannot directly minimize Eq. (10) with respect to n. Instead 
a scheme suggested by Kohn and Sham [10], is commonly employed for performing the minimization indirectly. This 
scheme starts by writing the minimization as 

5E ST S SV SU H SE XC ST 

As a consequence of Eq. (8), 5V/5n = u(r), the 'external' potential the electrons move in. The term SUn/Sn simply 
yields the Hartrcc potential J d 3 r' ^*J\ ■ The term v xc (r) = 5E xc /8n is universal exchange-correlation potential. 

Consider now a noninteracting system of particles moving in external potential v s (r). For this system the mini- 
mization condition is simply 

since there are no Hartree and xc terms in the absence of interactions. Comparing this with the previous equation 
we find that both minimizations have the same solution n s (r) = n(r) if v s is chosen to be 

v s (r) = v(r) + v H (r) + v xc (r). (13) 

Consequently, one can calculate the density of the interacting (many-body) system in potential v(r) by solving the 
equations of a noninteracting (single-body) system in potential v s (r) . In particular, the Schrodinger equation of this 
auxiliary system, 
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yields orbitals that reproduce the density n(r) of the original system 



&(r) = Ci &(r), (14) 



N 



n(r)^n s (r)=^|^(r)| 2 . (15) 

i 

Eqs. (13) to (15) are the celebrated Kohn-Sham equations. 

Since both vr and v xc depend on n, which depends on the fc, which in turn depend on v s , these equations should 
be solved till the self-consistency. 

From converged solution no, the total ground-state energy can be calculated using 
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E ° = f> - T j d * r j rfV n °| ( r r) _ n ° r / ( | ~ / <*Wr)n (r) + E xc [n \. 
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With the exact v xc (r) and i^n] all many-body effects are in principle included, which would be of tremendous 
importance for practical computational applications. 

Let us consider the following one-dimensional Schrodingcr equation 



& d2 ~\r ( ^ 



Mr) = eiMr), (17) 



If we have set of regularly ordered eigenvalues (e.g. = (i + ^)huj for harmonic oscillator), the local potential V(x) 
is given as a simple continuous function (V(x) = ^mujx 2 in this case). But let us consider that the infinite set 
of eigenvalues is chosen randomly. The mystery of how a one-dimensional integrable system can produce any given 
eigenvalues is resolved by adopting [11] the concept of a fractal potential V(x) which in the limit of infinite number 
of eigenvalues would have dimension higher than one. This principle can be illustrated on the example of nontrivial 
Riemann zeros, which were reproduced as an eigenvalues of Schrodingcr equation with potential having dimension 
d = 1.5 [12], or on infinite the sequence of prime numbers, where d — 1.8 [13]. It has been even suggested that in 
the limit of infinity number of eigenvalues the fractal dimension should converge to d = 2 [14]. In the Kohn-Sham 
equation exchange-correlation potential v xc (r) should be universal, which mean that the exact values of ground-state 
density and energy should be obtained for uncountable infinity of v(r) (all possible numbers of various nuclei at all 
possible positions in the space). Our suggestion is that in the analogy with the previous case this would be possible 
only when v xc (r) is a fractal function. 

it is possible to invert Kohn-Sham equation and obtain exchange-correlation potential for a given N-electron system 
from highly correlated electronic densities computed with traditional methods of quantum chemistry. As has been 
recently found, commonly used approximate functionals have potentials that are not even remotely close to the 
potentials reconstructed from very precise ab initio calculations on the two-electron atoms [15]. Nevertheless it should 
be stressed that such an exchange-correlation potential is constructed for a given particular system (atom) and hence 
is not free from an implicit dependence on the given particular external potential and obtained expression for v xc is 
exact only for given atom and for all other atoms or molecules constitutes only an approximation to the ground-state 
density. In this sense such an approach is analogous to the local-density approximation which similarly replaces exact 
universal exchange-correlation potential by that of homogeneous electron gas with the same local density. It is natural 
that potentials obtained in such a way are different from approximate functionals because these are constructed in 
order to be universal (although approximate) for all atoms and molecules. The exact universal potential would be 
also different from such a particular atomic potentials. 

For the reconstruction of genuine universal potential 10, 50 or maybe 100 atoms or molecules should be used 
simultaneously for iterative (till the Kohn-Sham equations would not reproduce exact densities for all atoms) re- 
construction of v xc (r) as done in the analysis of Riemann zeros potential (the simplest system for such a study are 
spherically symmetric atoms, a one-dimensional systems where v xc is just a function of r, the radial distance from 
nuclei). Nevertheless this is not an easy task because very high numerical precision is needed (in the reconstruction 
of the potential producing Riemann zeros or primes several hundreds of eigenvalues were needed to reveal fractal 
structure of a potential). Another possibility is to represents universal v xc (r) as a trial fractal function with several 
parameters to be fitted. 

Yet another possible way is to express the density-functional exchange and correlation potential in terms of the 
one-particle Green's function which in turn may be expressed in terms of the ground-state density Kohn-Sham orbitals 
[16]. In principle, this self-consistent loop yields the exact v xc (r). 

The concept of a fractal potential has been already successfully exploited in the study of protein conformational 
dynamics [17], where has been found that the fractal dimension of of the total potential energy as a function of time 
is almost independent of temperature, and increases with time, more slowly the larger the protein. Perhaps the most 
striking observation of this study is the apparent universality of the fractal dimension, which depends only weakly on 
the type of molecular system. 

It is interesting to note that both density functional theory [9,10] and Mandelbrot concept of fractals [18] emerged 
in the same time, about 40 years ago. We hope that our suggestion that fractals are also behind the density functional 
theory would be useful in searching for the exact universal exchange-correlation potential, the Holy Grail of many 
body physics. 
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